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LONG-TERM  GOAL 

The  long-term  goal  is  to  develop  an  efficient  methodology  for  determining  the  aerosol  optical 
properties  from  environmental  factors  (such  as  wind,  temperature,  humidity,  waves,  etc.)  in  the  lower 
part  of  the  marine  atmospheric  boundary  layer  (LP  MABL). 

OBJECTIVES 

Our  objective  is  to  develop  a  method  for  detennining  the  aerosol  microphysical  characteristics  in  LP 
MABL  from  the  backscattered  lidar  signal  (BLS).  This-year  objective  is  to  extend  the  previously 
constructed  inversion  algorithm  to  a  wider  class  of  aerosol  particle  size  distributions  (APSD),  and  to 
test  the  validity  of  the  inverse  method  by  experimental  lidar  data  obtained  by  the  group  of  Dr.  S. 
Sharma  of  the  Hawaii  Institute  of  Geophysics  and  Planetology  (HIGP). 

APPROACH 

Last  year,  we  developed  the  method  of  mean  ordinates  for  inverting  scarce  or/and  incomplete  aerosol 
optical  data  into  the  aerosol  microstructure.  The  method  proved  to  be  particularly  efficient  in  inverting 
the  aerosol  optical  characteristics  obtained  from  lidar  observations.  The  essence  of  the  method  is  as 
follows  [1]: 

•  The  APSD  is  described  by  a  sum  of  lognonnal  components  representing  small-sized,  medium¬ 
sized,  and  large  particles.  The  aerosol  parameter  intervals  are  intentionally  preset  so  wide  that  all 
known  types  of  atmospheric  aerosol  are  included  into  consideration. 

•  The  APSD  parameters  are  quantized  on  the  specified  intervals.  All  possible  combinations  of  the 
quantized  parameters  are  selected  by  exhaustion,  and  an  APSD  is  constructed  for  each 
combination. 

•  Optical  characteristics  of  each  constructed  model  are  calculated  and  compared  with  the  data  to  be 
inverted  (initial  data). 

•  The  models  whose  optical  characteristics  coincide  with  the  initial  data  within  preset  error  intervals 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  1999  2- REPORT  TYPE 

3.  DATES  COVERED 

00-00-1999  to  00-00-1999 

4.  TITLE  AND  SUBTITLE 

The  Determination  Of  The  Microphysical  Characteristics  Of  Aerosols  In 
The  Lower  Part  Of  The  Marine  Atmospheric  Boundary  Layer  From  The 
Structure  Of  The  Backscattered  Lidar  Signal 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROTECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Oregon  State  University, College  of  Oceanic  and  Atmospheric 

Sciences, 104  Oceanography  Admin  Bldg, Corvallis, OR, 97331 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  3S 

unclassified  unclassified  unclassified  Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


are  considered  as  acceptable  solutions. 

•  The  mean  ordinates  of  APSD  over  the  ensemble  of  all  acceptable  solutions  are  calculated. 

•  The  model  closest  to  the  curve  of  the  mean  ordinates  is  taken  as  the  most  probable  solution. 

The  method  was  tested  by  numerical  experiments  on  inverting  aerosol  optical  characteristics  obtained 
from  one-lidar  or  two-lidar  observations  [1].  The  results  proved  to  be  quite  satisfactory  for  smooth 
enough  APSDs.  In  this  respect,  as  well  as  in  terms  of  accuracy,  the  method  of  mean  ordinates  is 
comparable  to  the  constrained  inversion  method.  During  this  year,  we  have  studied  the  applicability  of 
the  method  to  a  wider  class  of  APSD  including  components  of  rather  narrow  distributions. 

The  method  of  mean  ordinates  has  been  employed  in  the  analysis  of  experimental  lidar  data  obtained  at 
HIGP  [2,3].  Horizontal  profiles  of  the  aerosol  spectral  attenuation  were  determined.  For  this  purpose, 
we  resorted  to  the  numerical  solution  of  the  basic  differential  lidar  equation  [4],  The  relationship 
between  the  aerosol  spectral  attenuation  c  and  backscattering  coefficient  (3  was  assumed  to  be  (3  = 
const  x  ok.  A  stable  solution  was  attained  by  specifying  the  boundary  condition  as  a  value  of  the 
spectral  attenuation  far  from  the  source.  It  was  calculated  by  the  slope  method  on  some  distance  range. 

APSD  was  determined  for  some  representative  points  of  a  profiles.  They  were  analyzed  in  terms  of  the 
behavior  of  particle  concentration  and  with  respect  to  the  shape  of  APSD  along  the  profiles. 

WORK  COMPLETED 

The  last-year  results  showed  that  the  method  of  mean  ordinates  works  quite  well  (with  an  error  of  15- 
25%)  for  one  or  two  lidars  provided  the  sought-for  APSD  is  smooth  enough.  The  smoothness  of  an 
APSD  curve  can  be  characterized  by  the  dispersion  of  APSD  components.  The  APSD  curve  is 
sufficiently  smooth  when  the  dispersion  is  greater  than  0.2.  When  the  dispersion  is  less  than  0.2 
(narrow  peaks),  data  of  one  or  two  lidars  are  inadequate  to  achieve  a  satisfactory  accuracy  of  the 
retrieved  APSD.  At  some  spectral  ranges,  the  retrieval  error  for  the  spectral  attenuation  may  be  as  high 
as  60%. 

We  applied  the  method  of  mean  ordinates  to  APSDs  containing  narrow  peaks  (0.1-0.15  dispersion). 

We  considered  three-lidar  data  at  0.55,  1.06,  and  10.591  pm,  and  four-lidar  data  at  0.55,  1.06,  4.7,  and 
10.591  pin.  It  turned  out  that  even  for  three  lidars,  the  inversion  results  are  promising:  the  height  and 
position  of  the  peaks  are  represented  quite  accurately.  The  results  on  the  peak  width  are  less  precise. 
The  retrieved  peaks  are  wider  than  the  initial  ones.  Notwithstanding  these  discrepancies  between  the 
initial  and  retrieved  APSD,  the  aerosol  optical  properties  are  represented  quite  accurately.  The  mean 
(over  the  whole  0.45-15  pin  spectral  range)  error  of  the  retrieved  spectral  attenuation  is  4%.  The 
maximum  error  does  not  exceed  16%  (at  7.2  pin). 

The  use  of  four  lidars  yields  still  better  results.  In  this  case,  the  difference  between  the  initial  and 
retrieved  spectral  attenuation  does  not  exceed  1%. 

The  results  of  the  inversion  of  three-lidar  data  obtained  by  the  method  of  mean  ordinates  are  shown  on 
Fig.la,b. 
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Fig.  1.  a):APSD;  1,  initial  APSD;  2,3,  minimal  and  maximum  ordinates;  4,  the  most  probable 
APSD;  b):  the  spectral  attenuation  o;  1,  initial  O;  2,  a  calculated  for  the  most  probable  APSD;  3, 

reference  values  of  O. 

We  calculated  profiles  of  the  spectral  attenuation  from  horizontal  lidar  data  at  0.532  pm.  The 
observations  were  carried  out  by  the  group  of  Dr.  Shanna  of  HIGP  on  the  Coconut  Island,  Hawaii,  in 
October  1997.  The  c  profiles  showed  significant  spatial  irregularities  of  dimensions  of  100-400  m.  In 
some  cases,  their  optical  thickness  was  as  high  as  0.27.  The  irregularities  took  place  against  the 
background  of  a  general  increase  of  a  in  the  seaward  direction. 

Using  the  method  of  mean  ordinates,  APSD  was  calculated  at  some  representative  points  on  the  c 
profiles.  Fig.  2  shows  a  profile  of  o  for  a  lidar  scan  on  October  16,  1997.  The  points  for  which  APSD 
were  detennined  are  marked  by  numbers. 
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Fig.  2.  Profile  of  the  spectral  attenuation  <7.  Six  points  with  error  bars  (±5%)  for  which  APSD  were 
retrieved  are  marked  by  numbers .  Circles  show  the  total  number  of  particles  (per  cm'3)  related  to  the 
numbered  points.  All  quantities  are  normalized  to  the  corresponding  values  at  point  2. 


In  order  to  unify  all  the  data  for  the  further  comparison,  they  were  reduced  to  those  at  the  point  2 
(corresponding  to  the  background  value  of  a).  Circles  mark  the  values  of  the  total  number  N  of 
particles.  It  is  seen  that  at  the  points  1,2,  and  6  variations  of  c  are  associated  with  variations  of  N 
while  the  shape  of  APSD  changes  only  slightly.  In  the  case  of  inhomogeneities  (points  3,4,5),  the 
increase  of  N  cannot  completely  account  for  the  increase  of  c.  Clearly  in  this  instance,  the  shape  of 
APSD  changes  significantly. 


Fig.  3.  Relative  deviations  of  APSD  at  the  numbered  points  from  APSD  values  at  point  2. 


Fig.  3  shows  relative  deviations  of  APSD  (in  reference  to  the  APSD  at  point  2)  normalized  to  one 
particle.  It  is  seen  that  at  the  points  3,4,5,  the  APSD  has  a  noticeable  hump  on  the  optically  active  0.5- 
10  pm  particle  size  interval,  as  compared  to  the  “background”  APSD  at  the  point  2.  This  accounts  for 
the  behavior  of  a  within  the  inhomogeneities.  In  contrast  to  this  picture,  deviations  of  APSD  at  the 
points  1  and  6  remain  within  25%  on  the  active  size  interval.  Noteworthy,  the  number  of  small 
particles  at  the  point  1  is  noticeably  higher  than  at  the  point  2.  In  other  words,  the  relative  content  of 
small  particles  recedes  from  the  coastline  in  the  seaward  direction. 

RESULTS 

1 .  It  was  established  that  the  use  of  greater  number  of  lidars  (more  than  two)  significantly  extends  the 
efficiency  of  the  method  of  mean  ordinates  for  inverting  lidar  data  into  APSD. 

2.  It  was  confirmed  that  significant  inhomogeneities  of  aerosol  optical  properties  take  place  in  the 
Hawaiian  coastal  zone.  Using  the  method  of  mean  ordinates,  we  estimated  the  magnitude  and 
character  of  these  inhomogeneities  in  terms  of  the  total  number  of  particles  as  well  as  with  regard 
to  variations  of  the  APSD  shape. 

IMPACT/APPLICATION 

The  method  of  mean  ordinates  proved  its  efficiency  in  inverting  lidar  data  into  the  aerosol  structure.  It 
has  been  tested  both  by  numerical  experiments  and  by  real  lidar  observation  data.  The  impact  of  these 
results  is  quite  significant.  Considering  the  fact  that  the  conventional  inversion  methods  fail  when 
applied  to  lidar  data  (because  of  a  limited  number  of  observation  wavelengths),  the  method  of  mean 
ordinates  is  currently  the  only  one  to  efficiently  utilize  lidar  information  for  marine  aerosol  studies. 

TRANSITIONS 

The  method  of  mean  ordinates  is  used  in  our  joint  work  with  the  experimental  group  of  Dr.  S.  Shanna 
of  the  Hawaii  Institute  of  Geophysics  and  Planetology  for  the  analysis  of  lidar  measurement  results  and 
for  studying  the  aerosol  structure  in  LP  MABL  in  the  Hawaiian  coastal  zone.  In  the  near  future,  the 
investigations  will  include  nephelometer  measurements  by  Dr.  A.  Hunt  of  the  Lawrence  Berkeley 
National  Laboratory.  The  method  of  mean  ordinates  has  been  employed  by  Dr.  P.  Bruscaglione, 
Instituto  Fisica  Superiore,  Italy,  for  studying  atmospheric  aerosol  near  the  sea  surface;  by  Dr.  O.V. 
Kopelevich,  Institute  of  Oceanology  of  the  Russian  Academy  of  Sciences,  for  hydrosol  studies;  by  Dr. 
J.  Lenoble,  Laboratory  of  Optics  of  the  Atmosphere,  Universite  des  Sciences  et  Technologie  de  Lille, 
France,  for  studying  atmospheric  aerosol. 

RELATED  PROJECTS 

The  approach  developed  in  the  method  of  mean  ordinates  for  inverting  scarce  or  incomplete  data  is 
being  used  in  our  NASA-funded  project  “The  improved  algorithm  for  estimating  the  atmospheric 
effect  in  space  measurements  of  the  chlorophyll  concentration.” 

REFERENCES 

1.  K.S.  Shifrin,  and  I.G.  Zolotov,  1999.  A  method  for  retrieving  the  aerosol  particle  size  distribution 
from  horizontal  lidar  data.  Submitted  to  Applied  Optics. 


2.  S.K.  Sharma,  B.R.  Lienert,  J.N.  Porter,  and  A.D.  Clarke,  1998.  Scanning  lidar  imaging  of  marine 
aerosol  fields  generated  by  breaking  waves.  In:  “Proc.  19th  International  Laser  Radar  Conf.  ”, 

pp. 673-676,  NASA  Langley  Research  Center,  Hampton  VA. 

3.  K.S.  Shifrin,  B.R.  Lienert,  J.N.  Porter,  and  S.K.  Sharma,  1999.  The  Determination  of 
Microphysical  Characteristics  of  the  Aerosol  in  the  Lower  Part  of  the  Marine  Atmospheric 
Boundary  Layer  from  the  Backscattered  Lidar  Signal.  Submitted  to  Remote  Sensing  of 
Environment. 

4.  J.D.  Klett,  1981.  Stable  analytical  inversion  solution  for  processing  lidar  returns.  Appl.  Optics, 
v.20,  211-220. 

PUBLICATIONS 

1.  K.S.  Shifrin,  1998.  Atmospheric  Aerosols,  Methods  of  Analysis.  In:  Encyclopedia  of 
Environmental  Analysis  and  Remediation,  v.  1,  446-505,  John  Wiley  &  Sons,  Inc.,  New  York. 

2.  I.M.  Levin,  K.S.  Shifrin,  1998.  Potential  for  Determining  the  Sea-Water  Absorption  Coefficient  By 
Satellite  Pulsed  Lidar.  Remote  Sensing  of  Environment,  v.  65,  1 105-11 11. 

3.  K.S.  Shifrin,  R.M.  Carlson,  and  I.G.  Zolotov,  1998.  Retrieval  of  the  Aerosol  Particle  Size 
Distribution  (APSD)  in  the  Lower  Part  of  the  Marine  Atmospheric  Boundary  Layer  (LP  MABL) 
from  Horizontal  Measurements  of  the  Atmospheric  Spectral  Transmittance.  Proc.  SPIE,  v.  3375, 
170-182. 

4.  A.M.  Kokorin,  K.S.  Shifrin,  1998.  Effect  of  the  Internal  Structure  of  Radially  Inhomogeneous 
Particles  of  a  Marine  Suspension  on  the  Scattering  Distribution  and  the  Degree  of  Plane 
Polarization  of  Scattered  Light.  Journal  of  Optical  Technology,  v.  65,  No  9,  26-32. 

5.  K.S.  Shifrin,  I.G.  Zolotov,  1999.  A  Method  for  Retrieving  the  Aerosol  Particle  Size  Distribution 
from  Horizontal  Lidar  Data.  Submitted  to  Applied  Optics. 

6.  K.S.  Shifrin,  B.R.  Lienert,  J.N  Porter,  and  S.K.  Shanna,  1999.  The  Detennination  of  Microphysical 
Characteristics  of  the  Aerosol  in  the  Lower  Part  of  the  Marine  Atmospheric  Boundary  Layer  from 
the  Backscattered  Lidar  Signal.  Submitted  to  Remote  Sensing  of  Environment. 

7.  K.S.  Shifrin,  1999.  The  angular  Structure  of  the  Upward  Radiance  under  the  Sea  Surface. 

Submitted  to  Journal  of  Geophysical  Research. 


